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Abstract 

The isoscalar giant dipole resonance structure in 208 Pb is calcu- 
lated in the framework of a fully consistent relativistic random phase 
approximation, based on effective mean-field Lagrangians with nonlin- 
ear meson self-interaction terms. The results are compared with recent 
experimental data and with calculations performed in the Hartree- 
Fock plus RPA framework. Two basic isoscalar dipole modes are iden- 
tified from the analysis of the velocity distributions. The discrepancy 
between the calculated strength distributions and current experimen- 
tal data is discussed, as well as the implications for the determination 
of the nuclear matter incompressibility. 
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The study of the isoscalar giant dipole resonance (IS GDR) might pro- 
vide important information on the nuclear matter compression modulus K nm . 
This, somewhat elusive, quantity defines basic properties of nuclei, super- 
novae explosions, neutron stars and heavy-ion collisions. The range of values 
of K nm has been deduced from the measured energies of the isoscalar giant 
monopole resonance (GMR) in spherical nuclei. The complete experimental 
data set on isoscalar GMR, however, does not limit the range of K nm to 
better than 200 — 300 MeV. Also microscopic calculations of GMR excita- 
tion energies have not really restricted the range of allowed values for the 
nuclear matter compression modulus. On one hand, modern non-relativistic 
Hartree-Fock plus random phase approximation (RPA) calculations, using 
both Skyrme and Gogny effective interactions, indicate that the value of 
K nm should be in the range 210-220 MeV ]I[ [2]. In relativistic mean-field 
models on the other hand, results of both time-dependent and constrained 
calculations suggest that empirical GMR energies are best reproduced by an 
effective force with K nm « 250 - 270 MeV || |, §. 

In principle, complementary information about the nuclear incompress- 
ibility, and therefore by extension about the nuclear matter compression mod- 
ulus, could be obtained from the other compression mode: giant isoscalar 
dipole oscillations. In first order the isoscalar dipole mode corresponds to 
spurious center-of-mass motion. The IS GDR is a second order effect, built 
on 3hu, or higher configurations. It can be visualized as a compression wave 
traveling back and forth through the nucleus along a definite direction: the 
"squeezing mode" [§, [7|. There are very few data on IS GDR in nuclei (the 
current experimental status has been reviewed in Ref. ||). In particular, 
recent results on IS GDR obtained by using inelastic scattering of a particles 
have been reported for 208 Pb §, and for 90 Zr, 116 Sn, 144 Sm, and 208 Pb Q. 
As in the case of giant monopole resonances, data on heavy spherical nuclei 
are particularly significant for the determination of the nuclear matter com- 
pression modulus: for example 208 Pb. However, recent experimental data 
on IS GDR excitation energies in this nucleus disagree: the centroid energy 
of the isoscalar dipole strength distribution is at 22.4 ± 0.5 MeV in Ref. M, 
while the value 19.3 ±0.3 MeV has been reported in Ref. pi. In the analysis 
of Ref. U , the " difference of spectra" technique was employed to separate the 
IS GDR from the high-energy octupole resonance (HEOR) in the 0° — > 2° 
a-scattering spectrum for 208 Pb. On the other hand, in the experiment on 
90 Zr, 116 Sn, 144 Sm, and 208 Pb of Ref. 0, the mixture of isoscalar L = 1 (IS 
GDR) and L = 3 (HEOR) multipole strength could not be separated by a 
peak fitting technique. Instead, the data were analyzed by a multipole anal- 
ysis of 1 MeV slices of the data over the giant resonance structure, obtained 
by removing the underlying continuum. 

In Ref. || it has been also pointed out that the experimental IS GDR 
centroid energies, and therefore the corresponding values of the nuclear in- 
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compressibility Ka, are not consistent with those derived from the mea- 
sured energies of the isoscalar GMR in 208 Pb. In the sum rule approach to 



the compression modes [ll|, two different models have been considered for 
the description of the collective motion: the hydrodynamical model and the 
generalized scaling model. The assumption of the scaling model leads to 
a difference of more than 40% between the values of the finite nucleus in- 
compressibility Ka, when extracted from the experimental energies of the IS 
GDR and the GMR in 208 Pb. A consistent value for K A can be derived from 
the experimental excitation energies, only if the two compression modes are 
described in the hydrodynamical model. The resulting value of Ka ~ 220 
MeV, however, is much too high, and in fact it corresponds to the nuclear 
matter compression modulus K nm , as derived from non-relativistic Hartree- 
Fock plus RPA calculations. This is not difficult to understand, since the 
expressions for Ka in both models were derived in the limit of large systems 
and consequently do not account for surface effects fll |. Both models, how- 



ever, are approximations to a full quantum description: the time-dependent 
Hartree-Fock or, equivalently, the RPA. Therefore, fully microscopic calcu- 
lations might be necessary, in order to resolve the apparent discrepancy be- 
tween the values of Ka extracted from the IS GDR and the GMR in 208 Pb. 

Non-relativistic self-consistent Hartree-Fock plus RPA calculations of dipole 
compression modes in nuclei were reported in the work of Van Giai and 



Sagawa fll2 |, and more recently in Refs. and [|Tj]. A number of different 
Skyrme parameterizations were used in these calculations, and the result is 
that all of them systematically overestimate the experimental values of the 
IS GDR centroid energies, not only for 208 Pb, but also for lighter nuclei. 
In particular, those interactions that reproduce the experimental excitation 
energies of the GMR (SGII and SKM*), predict centroid energies of the IS 
GDR in 208 Pb that are 4 — 5 MeV higher than those extracted from small 
angle a-scattering spectra. In Ref. fl4| effects that go beyond the mean-field 
approximation have been considered: the inclusion of the continuum and 
2p — 2h coupling. It has been shown that the coupling of RPA states to 
2p — 2h configurations, although it reproduces the total width, results in a 
downward shift of the resonance energy of less than 1 MeV with respect to 
the RPA value. It appears, therefore, that the presently available data on 
excitation energies of the compression modes in nuclei: the GMR and the IS 
GDR, cannot be consistently reproduced by theoretical models. 

In Ref. we have performed time- dependent and constrained relativis- 
ts mean-field calculations for the monopole giant resonances in a number of 
spherical closed shell nuclei, from 16 to 208 Pb. It has been shown that, in 
the framework of relativistic mean field theory, the nuclear matter compres- 
sion modulus K nm « 250 — 270 MeV is in reasonable agreement with the 
available data on spherical nuclei. This value is approximately 20% larger 
than the values deduced from non-relativistic density dependent Hartree- 
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Fock calculations with Skyrme or Gogny forces. In particular, among the 
presently available effective Lagrangian parameterizations, the NL3 effective 



force |ToJ with K nm = 271.8 MeV, provides the best description of the mass 
dependence of the GMR excitation energies. Preliminary calculations with 
the time-dependent relativistic mean-field model ||16|| , indicate that the NL3 
effective interaction, which reproduces exactly the excitation energy of the 
GMR in 208 Pb (14.1 MeV), overestimates the reported centroid energy of the 
IS GDR by at least 4 MeV. However, due to complications arising from the 
spurious center-of-mass motion, the time-dependent relativistic mean-field 
model computer code develops a numerical instability which prevents the 
precise determination of the IS GDR excitation energy. In the present analy- 
sis, therefore, we apply the relativistic random phase approximation (RRPA) 
to the description of the isoscalar dipole oscillations in 208 Pb. 

The RRPA represents the small amplitude limit of the time-dependent 
relativistic mean-field theory. Self-consistency will therefore ensure that the 
same correlations which define the ground-state properties, also determine 
the behavior of small deviations from the equilibrium. The same effective La- 
grangian generates the Dirac-Hartree single-particle spectrum and the resid- 
ual particle-hole interaction. Some of the earliest applications of the RRPA to 
finite nuclei include the description of low-lying negative parity excitations 
in 16 O ||17|| , and studies of isoscalar giant resonances in light and medium 



nuclei . These RRPA calculations, however, were based on the most sim- 
ple, linear a — u relativistic mean field model. It is well known that for a 
quantitative description of ground- and excited states in finite nuclei, den- 
sity dependent interactions have to be included in the effective Lagrangian 
through the meson non-linear self interaction terms. The RRPA response 
functions with nonlinear meson terms have been derived in Refs. [[TTJ pCfl , and 
applied in studies of isoscalar and isovector giant resonances. However, the 
calculated excitation energies did not reproduce the values obtained with the 
time-dependent relativistic mean-field model @, 16]. The reason was that the 



RRPA configuration spaces used in Refs. [19], |20[ did not include the negative 



energy Dirac states. In Ref. [21J it has been shown that an RRPA calculation, 
consistent with the mean-field model in the no — sea approximation, neces- 
sitates configuration spaces that include both particle-hole pairs and pairs 
formed from occupied states and negative-energy states. The contributions 
from configurations built from occupied positive-energy states and negative- 
energy states are essential for current conservation and the decoupling of 
the spurious state. In addition, configurations which include negative-energy 
states give an important contribution to the collectivity of excited states. 
In a recent study [22| we have shown that, in order to reproduce results of 
time-dependent relativistic mean-field calculations for giant resonances, the 
RRPA configuration space must contain negative-energy Dirac states, and 
the two-body matrix elements must include contributions from the spatial 
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components of the vector meson fields. The effects of the Dirac sea on the ex- 
citation energy of the giant monopole states have been also recently studied 
in an analytic way within the a — uj model ||23|| . 

In Fig. 1 we display the IS GDR strength distributions in 208 Pb: 

B T =\E1^^0 f ) = 1 -\(0 f \\Qj= \\l l )\ 2 , (1) 
where the isocalar dipole operator is 



Qlf = e E To (r 3 - w) YM, ifi), (2) 

and 



Hp 

i=l 



V = I < r 2 > ■ (3) 

The calculations have been performed within the framework of the self- 
consistent Dirac-Hartree plus relativistic RPA. The effective mean-field La- 
grangian contains nonlinear meson self-interaction terms, and the configura- 
tion space includes both particle-hole pairs, and pairs formed from hole states 
and negative-energy states. The choice of the dipole operator (0), with the 
parameter rj determined by the condition of translational invariance, ensures 
that the IS GDR strength distribution does not contain spurious components 
that correspond to the center-of-mass motion |i~2| . The strength distribu- 
tions in Fig. 1 have been calculated with the NL1 (K nm = 211.7 MeV) @], 
NL3 (1|] (K nm = 271.8 MeV), and NL-SH (K nm = 355.0 MeV) g| effective 
interactions. These three forces, in order of increasing values of the nuclear 
matter compressibility modulus, have been extensively used in the descrip- 
tion of a variety of properties of finite nuclei, not only those along the valley 
of /3-stability, but also of exotic nuclei close to the particle drip lines. In 
particular, in Ref. J|] it has been shown that the NL3 (K nm = 271.8 MeV) 
effective interaction provides the best description of experimental data on 
isoscalar giant monopole resonances. 

The calculated strength distributions are similar to those obtained within 
the non-relativistic Hartree-Fock plus RPA framework, using Skyrme effec- 
tive forces |T3, II4J. In disagreement with reported experimental results, all 



theoretical models predict a substantial amount of isoscalar dipole strength 
in the 8 — 14 MeV region. The centroid energies of the distributions in the 
high-energy region between 20 and 30 MeV, are 4 — 5 MeV higher than 
those extracted from the experimental spectra. It also appears that the cen- 
troid energies of the low-energy distribution do not depend on the nuclear 
matter incompressibility of the effective interactions. On the other hand, 
the IS GDR strength distributions in the low-energy region display the ex- 
pected mass dependence. We have also performed calculations for a number 
of lighter spherical nuclei, and verified that with increasing mass the cen- 
troid is indeed shifted to lower energy. When comparing with experimental 
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data, it should be pointed out that the usable excitation energy bite in the 
experiment reported in Ref. || was 14 — 29 MeV, and therefore a low-energy 
isoscalar dipole strength could not be observed. In this respect, somewhat 
more useful are the data from the experiment reported in Ref. fp70| , where 
spectra in the energy range m 4 < E x < 60 MeV have been observed. The 
results of an DWBA analysis of the experimental spectra, however, attribute 
the isoscalar strength in the 10 — 15 MeV region exclusively to the giant 
monopole (GMR) and giant quadrupole (GQR) resonances. It should be 
emphasized that a possible excitation of isoscalar dipole strength in this en- 
ergy region and its interference with the GQR cannot be excluded ||. 

In the high energy region the calculated dipole strength exhibits the ex- 
pected dependence on the nuclear matter compressibility modulus of the 
effective interactions (NL1, NL3, NL-SH). The centroid of the strength dis- 
tribution is shifted to higher energy with increasing values of K nm . These 
energies, however, are considerably higher than the corresponding experimen- 
tal IS GDR centroids || [TD[. Though, in order to precisely determine the IS 
GDR excitation energy from the experimental spectrum, the dipole strength 
has to be separated from the high-energy octupole resonance (HEOR), and 
this is not always possible QTO | . Using the NL3 effective interaction, we have 
calculated the octupole strength distribution. The centroid of the HEOR is 
found at « 22 MeV, well below the IS GDR main peak, but more than 2 
MeV above the experimental value for the HEOR centroid [[J. Incidentally, 
our calculated HEOR peak approximately coincides with the experimental 
value of the IS GDR centroid @. 

The IS GDR transition densities for 208 Pb are shown in Fig. 2. The tran- 
sition densities correspond to the NL3 strength distribution in Fig. 1. Since 
it appears that none of the effective interactions reproduces the experimental 
position of the IS GDR, the remainder of the present analysis will be only 
qualitative, and we choose to display only results obtained with the NL3 set 
of Lagrangian parameters. On the qualitative level, the other two effective 
interactions produce similar results. In Fig. 2. we plot proton (dot-dashed), 
neutron (dashed), and total (solid) transition densities for two representative 
peaks from Fig. 1: 10.35 MeV (a) is the central peak in the low-energy re- 
gion, and 26.01 MeV (b) is the energy of the main peak in the region above 
20 MeV. The transition densities for both peaks exhibit a radial dependence 
characteristic for the isoscalar dipole mode, and they can be compared with 
the corresponding transition densities in the scaling model, or with those 
which result from constrained calculations |ll|j. While for the high-energy 
peak the proton and neutron transition densities display an almost identical 
radial dependence, the pattern is more complicated for the peak at 10.35 
MeV. 

RPA calculations, therefore, predict the fragmentation of the isoscalar 
dipole strength distribution into two broad structures: one in the energy 
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window between 8 — 14 MeV, and the other in the high-energy region around 
ps 25 MeV. The position of the low-energy structure does not depend on 
the compressibility modulus, i.e. it does not correspond to a compression 
mode. Additional information on the underlying collective dynamics can 
be obtained through a study of transition currents. In Fig. 3 we plot the 
velocity fields for the two peaks at 10.35 MeV (a) and 26.01 MeV (b). The 
velocity distributions are derived from the corresponding transition densities, 
following the procedure described in Ref. p6| . The "squeezing" compression 



mode is identified from the flow pattern which corresponds to the high-energy 
peak at 26.01 MeV. The flow lines concentrate in the two "poles" on the 
symmetry axis at z ps ±2.5 fm. The velocity field corresponds to a density 
distribution which is being compressed in the lower half plane, and expands in 
the upper half plane. The centers of compression and expansion are located 
on the symmetry axis, at approximately half the distance between the center 
and the surface of the nucleus. It is obvious that the excitation energy of 
this mode will strongly depend on the compressibility modulus. The flow 
pattern for the lower peak at 10.35 MeV is very different. The flow lines 
describe a kind of toroidal motion, which is caused by the surface effect of 
the finite nucleus. The density wave travels through the nucleus along the 
symmetry axis. The reflection of the wave on the surface, however, induces 
radial components in the velocity field. Although it corresponds to dipole 
oscillations, this is not a compression mode. We have verified that also other 
dipole states in this energy region display similar velocity fields. 

In conclusion, the isoscalar giant dipole resonance in 208 Pb has been calcu- 
lated in the framework of the relativistic RPA, based on effective mean-field 
Lagrangians with meson self-interaction terms. The results have been com- 
pared with recent experimental data and with calculations performed in the 
Hartree-Fock plus RPA framework. While the results of the present RRPA 
study are consistent with previous theoretical analyses, they strongly disagree 
with reported experimental data on the position of the IS GDR centroid en- 
ergy in 208 Pb. This is a serious problem, not only because the disagreement 
between theory and experiment is an order of magnitude larger than for other 
giant resonances, but also because the present data on IS GDR are not con- 
sistent with the value of the nuclear incompressibility Ka derived from the 
measured excitation energy of the isoscalar GMR. This inconsistency could, 
perhaps, be explained by a possible excitation of isoscalar dipole strength in 
the low-energy window between 8 MeV and 14 MeV. Although predicted by 
all theoretical models, the low-lying IS GDR strength is not observed in the 
experimental spectra. From the analysis of the velocity fields, we have iden- 
tified two basic isoscalar dipole modes. The "squeezing" compression mode 
is found in the high-energy region at ps 26 MeV. The low-energy dipole mode 
does not correspond to a compression mode, and its dynamics is determined 
by surface effects. 
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Figure Captions 



• Fig.l IS GDR strength distributions in Pb calculated with the NL1 
(dashed), NL3 (solid), and NL-SH (dot-dashed) effective interactions. 

• Fig. 2 IS GDR transition densities for 208 Pb calculated with the 
NL3 parameter set. Proton (dot-dashed), neutron (dashed), and to- 
tal (solid) transition densities are displayed for the peaks at 10.35 MeV 
(a) and 26.01 MeV (b). 

• Fig. 3 Velocity distributions for the two isoscalar dipole modes in 
208 Pb calculated with the NL3 effective interaction. The velocity fields 
correspond to the two peaks at 10.35 MeV (a) and 26.01 MeV (b). 
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